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Abstract 

Ten commercial laboratories participated in an interlaboratory study to estab- 
lish the repeatability and reproducibility of compression strength tests conducted 
according to ASTM International Standard Test Method E9. The test employed a 
cylindrical aluminum AA2024-T351 test specimen. Participants measured elastic 
modulus and 0.2 % offset yield strength, YS{0.2 % offset), using an extensome- 
ter attached to the specimen. The repeatability and reproducibility of the yield 
strength measurement, expressed as coefficient of variations were cv r = 0.01 1 and 
cvr= 0.020 The reproducibility of the test across the laboratories was among the 
best that has been reported for uniaxial tests. The reported data indicated that us- 
ing diametrally opposed extensometers, instead of a single extensometer doubled 
the precision of the test method. Laboratories that did not lubricate the ends of 
the specimen measured yield stresses and elastic moduli that were smaller than 
those measured in laboratories that lubricated the specimen ends. A finite element 
analysis of the test specimen deformation for frictionless and perfect friction could 
not explain the discrepancy, however. The modulus measured from stress-strain 
data were reanalyzed using a technique that finds the optimal 6t range, and applies 
several quality checks to the data. The error in modulus measurements from stress- 
strain curves generally increased as the fit range decreased to less than 40 % of the 
stress range. 


Keywords: compression, testing, ILS, ASTM, E9, yield strength, elastic modulus 
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1 Introduction 


Compression testing is a conceptually simple method to establish the uniaxial stress- 
strain and mechanical behavior of materials. Because the test specimens can be right 
circular cylinders, they are easy to fabricate. The short gauge length of the test spec- 
imen sometimes makes it the only possible geometry for establishing uniaxial proper- 
ties, for example normal to the plane of a plate. Compression tests can also be used to 
establish the strength of brittle materials fTjEI that would be difficult to grip in tension. 
Kuhn 0 and Chait 0] have reviewed the methods for compression testing, and have 
demonstrated that although the test is conceptually simple, the user must overcome 
many experimental difficulties to translate the measured load and displacement curves 
into accurate stress-strain behavior. 

Much of the literature, see Table [I] on implementing the compression test focuses 
on tests to large strains ||5]|7l[8l flOl [T2lTl3l 11711181 . to establish work hardening behavior 
beyond the strains at which tension test specimens neck. Other studies are devoted 
to modelling processes such as upsetting 

Emin] mum. A common theme in 
these and other studies has been to quantify the effect of friction USEE] or optimal 
lubricant 0 0 nu m and barreling SHU El El on interpreting the stress-strain 
behavior of the test specimen. 

ASTM International standard test method E9, ill first established in 1924, is a con- 
sensus standard for conducting compression tests to establish the strength of materials. 
It contains methods for testing cylindrical test specimens as well as methods for test- 
ing sheets with lateral support. The standard contains requirements for calibrating and 
qualifying the testing machine including extensometers, and aligning the fixtures and 
test specimen. It recommends, but does not require, specific test specimen geometries 
and lubricants. In addition, it suggests several fixture designs. It does not require the 
use of an extensometer in contact with the test specimen to measure strain, and most 
studies that reference E9 infer the test specimen strain and stress from the displacement 
of the actuator. 

Although the first version of E9 was released eight decades ago, its precision has 
never been formally evaluated, as required m. Such evaluations require a formal 
interlaboratory study, which can benefit both end user and testing laboratories. Users 
need the results of interlaboratory tests to determine the uncertainty that should be 
associated with the value of a material parameter obtained using a test method. Labo- 
ratories that employ test methods use interlaboratory studies to identify the deficiencies 
in their test methods and improve their implementation of them. 

This manuscript reports the results of an interlaboratory study to establish the preci- 
sion of ASTM International standard test method E9 for determining the yield strength 
and elastic modulus in compression. The results of this study were incorporated into 
the Precision and Bias statement of E9 in 2009. This report goes beyond the research 
report li20l that documents the calculation of the precision by comparing the results to 


other interlaboratory studies of uniaxial test methods. Sec. 4.1 analyzing some of the 


possible sources for the variability,Sec. 4.2 and presenting a method to evaluate the 
elastic modulus measured in the test, Sec|4.3| 
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Table 1 : Summary of literature that analyzes the compression test. 


Reference 

Content 

^max 

Materials 

Lubricant 

Banerjee, 19850 

b 

1.6 

A1 

teflon, M 0 S 2 , oil. 





none 

Carter, 19850 

F 

0.35 

A1 

M 0 S 2 , none 

Chait, 19750 

r 

NA 



Cook, 19450 

f 

0.7 

Cu 

none 

Gunasekera, 1982(51 

f 

1.2 

1022 steel 

teflon 

Hsii, 19690 

bfF 

NA 

Cu 

teflon 

Kamaluddin, 2007 (TO) 

F 

0.8 

A1 

grease 

Kobayashi, 19701m 

b 

1.4 

1040 steel 

graphite 

Lovato, 1992H1 

f 

1 

Al, Nb, brass, 

M 0 S 2 , teflon. 




steel 

BN, none 

Male, 1966(131 

f 

0.8 

Al, Ti, brass 

graphite, lanolin, 
paraffin 

Mescall, 1983(151 

F 

NA 

4340 steel 


Papirno, 1 983 lfT5l 

0 

NA 

steel 

teflon 

Ray, 1983(161 

f 

NA 

steel 

teflon 

Schey, 1982(171 

b 

1 

1020 steel, 6061 
Al 

M 0 S 2 , teflon 

Woodward, 1977U51 

f 

1.2 

steel 

teflon 


e m ax • maximum strain in test 

NA : not available 

Key to content of reference 

b: barreling analysis 

f: friction analysis 

F: finite element analysis 

r : review article 


2 Experimental Procedure 

This interlaboratory study followed the methods of ASTM E691 BTI . and uses statis- 
tical terms in accord with ASTM E177 J22). 

2.1 Participants 

Using the ASTM International (23l and American Association for Laboratory Accred- 
itation (24l laboratory directories, the organizers contacted and discussed the interlab- 
oratory study with twenty-five possible laboratory participants. From this original list 
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Table 2: List of participants in this study. 


Participant 

URL 

AADFW, Inc, Euless, Tx 

http : / / www . aadfwinc . com 

Alcoa, Pittsburgh, Pa 

http : / / www . alcoa . com 

Exova, Glendale Hts, IL 

http : / /www . exova . com 

Dickson Testing Company Inc., 
South Gate, Ca 

http : / / www . dicks onte sting . com/ 

Imperial College Mechanical 
Engr., London, England 

http : / / www . imperial . ac . uk/ 

MAR-TEST, Inc. (Cincinnati), 
Cincinnati, Oh 

http : / / www . mar- test . com/ 

MAR-TEST, Inc. (Stuart), Stuart, 
FI 

http : / / www . mar- test . com/ 

Metcut Research Inc., Cincinnati, 
Oh 

http : / / www . metcut . com 

Stork Climax Research Services, 
Wixom, MI 

http : / / www . storksmt . com/ cr s 

Westmoreland Mechanical 
Testing, Youngstown, Pa 

http : / / www . wmtr . com 


thirteen laboratories agreed and were able to participate, and ten ultimately completed 
the test program. Table [2] 

2.2 Instructions and method 

The participants followed ASTM Standard Test Method E9 0 to establish the elastic 
modulus, E, and 0.2 % offset yield strength, FS(0.2 % offset). At the time of the study, 
the version of E9 in use was E9-89a, but the only non-editorial difference between E9- 
89a and the current version, E9-09 was the addition of the precision statement to the lat- 
ter, which was the purpose of the interlaboratory study. The participants also returned 
electronic traces of the stress-strain curves to the organizers. All compression fixtures 
were required to be qualified according to ASTM E9 m Section 6.6 using at least five 
of supplied test specimens, unless the participant had already qualified the test setup ac- 
cording to Section 6.6. The participants conducted the compression tests using at least 
one extensometer at a nominal strain rate of de/dt = 0.005 min -1 = 8.33 x 10 -5 s -1 . 
No laboratory reported strain measured from actuator displacement. Each participant 
reported ten items as required in sections 10.1.1-10.1.9 and 10.1.13 of standard method 
E9: Material (test specimen ID), configuration description, test specimen dimensions 
as tested, fixture and lubricant description, testing machine description, speed of test- 
ing (required in section 8.7; report actual value), stress-strain diagram, modulus of 
elasticity, E, yield strength, YS( 0.2 % offset), and any anomalies. 
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Figure 1: Schematic of the location of test specimen blanks in the original plate. 


2.3 Test Specimen 

The material tested was aluminum alloy AA2024-T351, which is solution heat-treated 
and stress-relieved by controlled stretching. It was supplied as a plate with thickness, 
t = 22.2 mm (0.875 in), from which the organizers cut and distributed test specimen 
blanks, which were distributed throughout the plate. Figure[T]describes the location of 
the test specimens in the original plate and their numbering scheme. Only test speci- 
mens in rows 1 and 2 were used in this study. Test specimens number 1-50 came from 
row 1; test specimens 51-100 came from row 2. Figure[2]shows the test specimen with 
dimensions and tolerances. The drawing provided to the participants showed the di- 
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Figure 2: Dimensions of the test specimen used. 


mensions in non-SI units. Participants machined their own test specimens from sawed 
blanks that the organizers supplied. 

Test specimens were tested with the loading axis transverse to the rolling direc- 
tion in the plane of the plate, the so-called long-transverse (“LT”) orientation, Mil- 
Handbook-5J li25l Figure 1.4.12.3(a). The test specimen ID takes the form “LT-NN-L- 
X” where “FT” identifies the orientation relative to the rolling direction, “NN” identi- 
fies the test specimen number (1-250), “L” identifies the laboratory that received the 
test specimen blank (A-K) and “X” identifies the test order of the test specimen at a 
given laboratory. 

2.4 Compression fixtures 

The participants used a variety of loading fixtures, some of which involved a sub-press 
mounted in the testing machine to improve alignment. Others used adjustable platens 
that were aligned and then locked. Three laboratories used diametrically opposed ex- 
tensometers instead of a single extensometer. No laboratory reported strain from ac- 
tuator displacement. Table [3] describes the specimen fixturing, alignment capability, 
number of extensometers, and the lubricant used in each laboratory. 
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Table 3: Descriptions of the compression setup, extensometers and gauge lengths, G, 
and lubricants for all laboratories. 


Lab 

test setup 

Extensometer 

Lubricant 

A 

No sub-press 

Two, opposed, 
class B-l 
G = 25.4 mm 

Not reported 

B 

Sub-press 

Single, class not 
reported 
G = 12.7 mm 

Molybdenum disulfide 

C 

No sub-press, no 

Single, class 

None 


spherical seat 

not-reported 
G = 25.4 mm 


E 

sub-press 

Single, class B-l 
G = 25.4 mm 

Teflon tape 

H 

No sub-press 

Single class B-2, 
G = 12.7 mm 

Molybdenum disulfide 

I 

No sub-press 

Two opposed, 
class B-2, 

G = 25.4 mm 

Not reported 

J 

Precision ground 

Two opposed. 

None 


sub-press, aligned to 

class B-l 



closer tolerance than 
required for specimen, no 
spherical seat 

G = 25.4 mm 


K 

No sub-press. Platens 

Single, class 

Molybdenum disulfide 


aligned, shimmed and 

B-2, 



then locked. No spherical 

G = 12.7 mm 



seat 



L 

Sub-press 

Single, class B-2 
G = 25.4 mm 

WD-40 

M 

Compression platens 

Single, class 

None 


mounted in aligned 

B-2, G not 



hydraulic grips 

reported 



3 Results 

3.1 Stress-strain behavior 

Figure [3] plots the engineering stress-strain curves by laboratory. For convenience, 
compressive strains and stresses are plotted as positive, rather than negative, values. 
In each case the stress-strain curve was shifted along the strain axis so that a linear 
fit to the stress-strain ( S — e) data in the range (50 < S < 175) MPa intercepts 
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Figure 3: Engineering stress-strain curves. 


the origin. Figure [4]presents four views of the same stress-strain curves over different 
strain ranges. 

Table [4] summarizes the reported modulus, E, and yield strength, YS(0.2 % offset), 
for the laboratories. Figure [5] plots the reported elastic modulus, E, by laboratory. 
Figure [6] plots the reported 0.2 % offset yield strength, YS(0.2 % offset), by laboratory. 


3.2 Testing rates 

Figure[7]plots the strain as a function of time for the laboratories that reported time data. 
The interlaboratory instructions did not specify a control mode for the test. Four lab- 
oratories (C,H,K,F) conducted the test in strain control from the extensometer signal, 
as indicated from the constant slope of the strain-time plot. Figure [7] Three labora- 
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Figure 4: Engineering stress-strain curves: (a) complete behavior, (b) low-strain be- 
havior, (c) yield behavior, and (d) plastic behavior. Dashed lines show the accepted 
value for the modulus of 2024-T351, E = 75.2 GPa, the 0.2 % offset yield strength 
determination, and the e = 0.008 total elongation. 


tories (A,E,J) conducted the test in position control, as indicated from the changing 
slope in the strain-time plot. Figure [8] plots the strain rates in the elastic and plas- 
tic portions of the stress-strain curves calculated by linear regression. The plastic 
strain rate was calculated for strains greater than the strain at the 0.2% offset yield 
strength: eysoo 2 < e < e max - The elastic strain rate was calculated in the range 
0 < e < (eysoo 2 — 0.002). Laboratory J conducted the test in position control and 
set the elastic strain rate close to the specified value. 
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Table 4: Reported data 


Lab 

Sample 

E 

GPa 

YS 

MPa 

Lab 

Sample 

E 

GPa 

YS 

MPa 

Lab 

Sample 

E 

GPa 

YS 

MPa 

A 

LT-85 

75.6 

339.0 

H 

LT-10 

77.2 

345.0 

K 

LT-24 

75.8 

343.0 

A 

LT-1 

72.9 

340.0 

H 

LT-44 

73.8 

347.0 

K 

LT-42 

80.7 

348.0 

A 

LT-40 

76.8 

341.0 

H 

LT-55 

73.1 

347.0 

K 

LT-51 

87.6 

351.0 

A 

LT-26 

75.1 

342.0 

H 

LT-63 

77.2 

347.0 

K 

LT-56 

84.1 

351.0 

A 

LT-66 

73.8 

342.0 

H 

LT-76 

75.8 

347.0 

K 

LT-5 

80.0 

352.0 

A 

LT-78 

74.7 

342.0 

H 

LT-91 

77.2 

347.0 

K 

LT-1 3 

76.5 

354.0 

A 

LT-82 

76.3 

343.0 

H 

LT-1 9 

78.6 

354.0 

K 

LT-8 3 

NA 

NA 

B 

LT-28 

73.1 

348.0 

I 

LT-43 

75.0 

348.0 

L 

LT-79 

73.0 

348.0 

B 

LT-39 

76.5 

349.0 

I 

LT-57 

76.0 

349.0 

L 

LT-45 

84.2 

349.0 

B 

LT-54 

75.0 

349.0 

I 

LT-65 

75.0 

349.0 

L 

LT-69 

78.4 

352.0 

B 

LT-9 

75.0 

353.0 

I 

LT-21 

76.0 

350.0 

L 

LT-1 2 

84.2 

355.0 

B 

LT-80 

83.4 

356.0 

I 

LT-34 

75.0 

350.0 

L 

LT-36 

75.7 

356.0 

C 

LT-68 

72.0 

318.0 

I 

LT-74 

75.0 

350.0 

L 

LT-6 

77.5 

356.0 

C 

LT-84 

72.8 

335.0 

I 

LT-70 

76.0 

351.0 

L 

LT-92 

85.8 

363.0 

C 

LT-90 

71.5 

336.0 

J 

LT-4 

74.5 

340.5 

M 

LT-1 8 

69.0 

338.0 

C 

LT-1 6 

72.4 

338.0 

J 

LT-89 

74.2 

341.6 

M 

LT-1 1 

72.0 

342.0 

C 

LT-77 

71.4 

338.0 

J 

LT-94 

74.4 

341.6 

M 

LT-59 

71.0 

342.0 

C 

LT-27 

70.7 

339.0 

J 

LT-73 

74.5 

342.1 

M 

LT-62 

70.0 

342.0 

C 

LT-58 

72.9 

341.0 

J 

LT-37 

74.8 

342.5 

M 

LT-8 6 

72.0 

342.0 

E 

LT-25 

68.3 

342.2 

J 

LT-23 

74.8 

343.1 

M 

LT-33 

72.0 

343.0 

E 

LT-31 

67.0 

346.8 

J 

LT-1 4 

74.1 

343.4 

M 

LT-81 

68.0 

345.0 

E 

LT-87 

69.3 

347.4 









E 

LT-71 

68.6 

349.0 









E 

LT-1 7 

70.5 

351.1 









E 

LT-8 

70.9 

351.4 









E 

LT-41 

70.6 

355.8 










Notes: 

YS=YS( 0.2 % offset) 


4 Analysis and discussion 

4.1 Expected variability 

Table [5] summarizes the mean and standard deviation, sd, and coefficient of variation, 
cv, of the reported n measurements of the elastic moduli and yield strengths. Table [6] 
defines the statistical parameters used in this section. 

The results of an interlaboratory study contain variability that arises within a given 
laboratory and variability that arises between laboratories. The terms repeatability and 
the reproducibility, denoted by subscripted “r” and “R" are often used to differenti- 
ate between the sources J22J. In a general sense, repeatability characterizes the ability 
of an individual laboratory to repeat measurements, while reproducibility characterizes 
the ability of an individual laboratory to achieve the global mean or accepted value. For 
example, if the results from an individual laboratory are tightly grouped, their repeata- 
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Figure 5: Reported elastic modulus, E. Dashed lines show the average reported mod- 
ulus and a commonly accepted value of the modulus ( E=15.2 GPa) from Mil-HDBK- 
5J l25l and the average of all the measurements, E = 74.9 GPa. 


bility is high, but their mean value may still deviate significantly from the accepted 
value, in which case their reproducibility is low. 

Both the variability in the material and the variability of the test method influ- 
ence the repeatability. The excellent repeatability of the measurements for both E and 
75(0.2 % offset) in laboratories A, I, and 1 , Table [5] shows that the material variability 
in this study is quite low, and therefore the implementation of the test method by the 
individual laboratories is the major contributor to repeatability. 

Figure [9] and Table [7] compare the results for repeatability and reproducibility for 
this study to the results of other interlaboratory studies of tensile tests at room 126]|27]| 
and elevated temperatures (28] to establish the 0.2 % offset yield strength, T5(0.2 % 
offset). In Figure[9]the within-laboratory and between-laboratory standard deviations. 
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Figure 6: Reported 0.2 % offset yield strength, YS( 0.2 % offset). The dashed line is the 
average of the results from all laboratories. 


s r , Eq. [4] and Sr, Eq. [7] see Table [6] are divided by the mean value and expressed as 
their respective coefficients of variation, cv r , Eq.[5] and cvr, Eq. [8] 

The repeatability coefficient of variation of the yield strength measurements in this 
study, derived from the average of the standard deviations of the individual laborato- 
ries, is cv r = 0.0111, Eq. [5] while the coefficient of variation of the reproducibility, 
cvfj = 0.0197, Eq. [8] is about twice as large. The repeatability and reproducibility 
of compression tests established in this study are among the best measurements of all 
reported uniaxial measurements l26l |27I [28] | . 
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Figure 7: Reported strain as a function of time. 


4.2 Effect of test method implementation 

Because the laboratories were free to implement the test method in different ways, 
while still complying with the standard, their result can be used to learn about the effect 
of different aspects of test method on the measurement. In particular, the results can be 
examined to determine the effects of the number of extensometers and lubrication on 
the quality of the measurement. 

4.2.1 Effect of number of extensometers 

Some laboratories used a single extensometer, while others used two opposed exten- 
someters to measure the strain. Figure [lOjplots the reported modulus and yield strength 
for the laboratories, grouped by the number of extensometers used in the determination. 
The three laboratories. A, I, J, with the best repeatability and reproducibility of E and 
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Figure 8: Calculated (a) elastic and (b) plastic strain rates for laboratories that re- 
ported time data. Symbols differentiate between position-control and strain-control 
tests. Dashed line indicates the requested rate, de/dt = 8.33 x 10 -5 1/s. 


Monday 27 th September, 2010: WERB Version: Technical Note 


15 


0.10 - 
0.08 - 

O 


alloy steel O 

aluminum A 

stainless steel <0> 
steel 

superalloy □ 

This Study 

0.06 - 

o 



- 

0.04 - 




- 

0.02 - 

A 

° 0 
o 

□ 

o 

o 

- 

0.00 - 




- 


200 

400 

600 800 

1000 



Yield Strength/ MPa 




(a) Repeatability 


0.10 - 
0.08 - 

o 

o 


alloy steel O 

aluminum A 

stainless steel O 
steel 

superalloy □ 

This Study 



O 



0.06 - 




- 

0.04 - 

c 

o o 

O o 

o 

o 

- 

0.02 - 

s 

A 

A 

A 


O 

0.00 - 




- 


200 

400 

600 800 

1000 


Yield Strength/ MPa 


(b) Reproducibility 


Figure 9: Literature data ll26l [281 [27] on (a) repeatability (within-laboratory) and (b) 
reproducibility (between-laboratory) of yield strength, KSY0.2 % offset). Table[7]sum- 
marizes the data. 
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Table 5: Summary statistics for the modulus and yield strength data. 


Lab 

n 

E sd (E) 
GPa GPa 

cv(E) 

75 

MPa 

sd(75) 

MPa 

cv(YS) 

A 

7 

75.0 

1.4 

0.018 

341.3 

1.4 

0.004 

B 

5 

76.6 

4.0 

0.052 

351.0 

3.4 

0.010 

C 

7 

72.0 

0.8 

0.011 

335.0 

7.7 

0.023 

E 

7 

69.3 

1.4 

0.021 

349.1 

4.3 

0.012 

H 

7 

76.1 

2.0 

0.026 

347.7 

2.9 

0.008 

I 

7 

75.4 

0.5 

0.007 

349.6 

1.0 

0.003 

J 

7 

74.5 

0.3 

0.004 

342.1 

1.0 

0.003 

K 

7 

80.8 

4.5 

0.056 

349.8 

3.9 

0.011 

L 

7 

79.8 

4.9 

0.062 

354.1 

5.1 

0.014 

M 

7 

70.6 

1.6 

0.023 

342.0 

2.1 

0.006 


Note: YS= YS( 0.2 % offset) 


Parameter 

Grand Average 

s r 

cv r 

sr 

cv R 

E 

T5(0.2 % offset) 

E= 75.0 GPa 

75(0.2 % offset)= 346.2 MPa 

2.7 GPa 

3.8 MPa 

0.036 

0.011 

4.4 GPa 
6.8 MPa 

0.059 

0.020 


Table[6]defines the parameters. Because laboratory B tested only five specimens instead of s 


the mean values shown in Figures p 
the average of cell averages, Table [6 


6] |10[ and|12|differ from the Grand Average, also known as 
EqUl 


YS( 0.2 % offset) are the only ones that used a system of two-opposed extensometers in- 
stead of a single extensometer. In addition, the measured elastic modulus in these three 
laboratories was the closest to the accepted f25l value E = 75.2 GPa, Figure[lOji. Note 
that the repeatability of the modulus and yield strength from the labs that used a single 
extensometer is worse (i.e. larger variability), but also that for an individual laboratory, 
the values usually all lie above or below the grand average. Using two extensometers 
will tend to average the effect of non-axial loading, since during bending, one side of 
the test specimen will be displaced more, and the other less. Because the values for 
an individual, single -extensometer laboratory are usually displaced to one side of the 
global mean value points toward the alignment of the compression fixture, rather than 
machining of the test specimen as the source of the non-axiality. If the ends of the test 
specimen were not parallel, the measured values would tend to encompass the mean 
value, because no special relationship exists between the test specimen and the loading 
fixture. Conversely, the extensometer is usually mounted on the test specimen with a 
fixed relation to the orientation of loading fixture, so bending induced by a slightly mis- 
aligned loading fixture is always in the same geometric relation to the extensometer, 
which will tend to always over- or under-estimate the strain. 


Monday 27 th September, 2010: WERB Version: Technical Note 


17 


Table 6: Definitions of repeatability and reproducibility statistics. 


Formula 

n 

P 

Y 

1 


Description 

number of tests in a laboratory, n ss 7 
number of laboratories, p = 10 
individual test result in lab j 

1 < i < n ; 1 < j < P 


Yj 

f 



average test result in lab j (cell average) (1) 
average of cell averages (grand average) (2) 


A 


i n 


yu 


E (^) 2 


standard deviation measured in lab j (3) 
repeatability (within-lab) standard deviation (4) 


CV r =^et 
Y 

coefficient of variation within laboratories 

(5) 


pljEm r > 2 

standard deviation of cell averages 

(6) 

SR =v 

^4 + s '( n ) 

reproducibility standard deviation ((between-lab) 

(7) 

Sr 

CVr = — 

coefficient of variation between laboratories 

(8) 


Y 


4.2.2 Effect of Lubricant 

The standard allows the laboratory to choose whether to lubricate the ends of the test 
specimen during the test. Figure [IT] shows the engineering stress-strain curves, identi- 
fied by the use or omission of lubricants. Several laboratories did not report lubricant 
use. Table [3] identifies the lubricants used. Three laboratories used no lubricant, and 
their stress-strain curves are the lowest in the collection. Figure [12] plots the reported 
modulus, E, and yield strength, KSY0.2 % offset), identified by lubricant use. The aver- 
age yield strength from the three laboratories that did not use lubricant was 10.6 MPa 
lower (3.1 %) than the average of the laboratories that lubricated the specimen. Simi- 
larly the average reported elastic modulus, E, from the three laboratories was 4. 1 GPa 
(5.7 %) smaller. Table[8] an analysis of variance, shows that the reported yield strength, 
KSY0.2 % offset), and modulus, /Y,are both significantly different for the lubricated and 
unlubricated cases. 
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Table 7: Literature data 11281 [26 , 27 1 for repeatability and reproducibility in mechanical 
testing. 


Source 

Material 

T 

°C 

YS 

MPa 

cv r 

CVR 

[261 

EC-H19 aluminum 

20.0 

158.4 

0.0210 

0.0210 

El 

2024-T351 aluminum 

20.0 

362.9 

0.0140 

0.0150 

Cl 

AISI A 105 steel 

20.0 

402.4 

0.0140 

0.0250 

Ea 

SS316 stainless steel 

20.0 

480.1 

0.0140 

0.0410 

Mi 

Inconel 600 

20.0 

268.3 

0.0090 

0.0220 

ESI 

SAE 51410 steel 

20.0 

967.5 

0.0090 

0.0160 

E3 

AA5754 A1 

20.0 

105.7 

NA 

0.0160 

E3 

AA5 1802-0 

20.0 

126.4 

NA 

0.0100 

E3 

AA6016-T4 

20.0 

127.2 

NA 

0.0110 

E3 

DX56, low-carbon steel 

20.0 

162.0 

NA 

0.0230 

Ell 

HR3 steel plate 

20.0 

228.6 

NA 

0.0410 

E3 

ZStE 180 steel 

20.0 

267.1 

NA 

0.0500 

E2 

S355 steel plate 

20.0 

427.6 

NA 

0.0310 

E3 

SS316L stainless steel 

20.0 

230.7 

NA 

0.0350 

E3 

X2CrNil8-10 stainless steel 

20.0 

303.8 

NA 

0.0330 

El 

X2CrNiMol8-10 

20.0 

353.3 

NA 

0.0390 

E3 

30NiCrMol6 high strength steel 

20.0 

1039.9 

NA 

0.0100 

E3 

Nimonic 75 CRM 661 

20.0 

302.1 

NA 

0.0180 

E3 

Nimonic 75 CRM 661 

20.0 

298.1 

NA 

0.0200 

El 

SS304 stainless steel 

316.0 

127.3 

0.0580 

0.1040 

El 

Low-carbon steel 

316.0 

236.4 

0.0110 

0.0350 

El 

2.25 Cr IMo steel 

316.0 

454.7 

0.0230 

0.0330 

El 

A286 stainless steel 

316.0 

699.5 

0.0280 

0.0300 

El 

SS304 stainless steel 

593.0 

101.4 

0.0770 

0.0990 

El 

Low-carbon steel 

593.0 

133.0 

0.0390 

0.1040 

El 

2.2.5Cr IMo ferritic steel 

593.0 

337.9 

0.0260 

0.0670 

El 

A286 stainless steel 

593.0 

642.2 

0.0260 

0.0350 

This study 

A1 2024-T351 

20.0 

346.2 

0.0111 

0.0197 


Notes: 

YS= YS{0.2 % offset) 
NA = data not available 


Since most laboratories used an extensometer with a gauge length G = 25.4 mm on 
the test specimen. Fig. [2] with length l « 48 mm and diameter d = 12.7 mm, the end 
of the test specimen is approximately one test specimen diameter away from the exten- 
someter contact point. That distance complies with the requirements of E9, and should 
minimize the effect of the non-uniform stress caused by friction at the specimen ends. 
No study has examined the sensitivity of the strain measurement to the non-uniform 
stress at the end of the specimen, but the existence of non-uniform stresses caused by 
friction is well known l29l . Finite element analysis (FEA), using the commercial finite- 
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Table 8: Analysis of variance tables for assessing the effect of lubrication on modulus 
and yield strength. 


% Table generated on 2010-09-22 19:05:27 
Analysis of Variance Table 

Model 1: E " 1 

Model 2 : E Lubricated 

Res.Df RSS Df Sum of Sq F Pr(>F) 

1 52 1180.16 

2 51 971.31 1 208.85 10.966 0.00171 ** 

Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1 

Analysis of Variance Table 

Model 1: YS02 “ 1 

Model 2 : YS02 ~ Lubricated 

Res.Df RSS Df Sum of Sq F Pr(>F) 

1 52 2655.5 

2 51 1223.9 1 1431.6 59.653 3.939e-10 *** 

Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1 


Key to parameters in table: 

E: E, elastic modulus 

YS02: YS(0.2 % offset), 0.2 % offset yield strength 

Data from laboratories that did not report the lubrication are omitted. Model 1 is the 
mean value of of the parameter, and model 2 computes the mean value for the lubricated 
(“Yes”) and unlubricated (“No”) cases. 


element modeling software Abaqus/Standard 6.10155). was used to analyze the effect 
of friction and gauge length on the response of the specimen. By exploiting the axial 
and radial symmetry of the specimen and platens, it was only necessary to model the 
upper-right quadrant the specimen with radius r = 6.3 mm and length l = 24.25 mm. 
The specimen mesh. Figure [13] was constructed with 1542 quadrilateral four-node 
axi-symmetric (CAX4) elements. The platens were modeled as non-deforming, rigid 
planes of infinite elastic modulus. Roller boundary conditions were applied on the axial 
centerline and the specimen midplane. The specimen material was modeled as a ho- 
mogeneous elastic -plastic time-independent material exhibiting strain hardening. The 
plastic deformation was modeled by the .l-> (or von Mises) flow plasticity theory with 
isotropic hardening. 130) The uniaxial stress-strain curve was taken from specimen LT- 
14 of the experimental results, down-sampled to 58 points, with the Young’s modulus 
set to the accepted value E = 75.2 GPa, and Poisson’s ratio=// = 0.33. This material 
model yields at a true stress a = 190.4 MPa at a true strain e = 0.0026, and not at the 
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value of the 0.2 % offset value 75(0.2 % offset)=343.4 MPa. All of the stress-strain 
curves, Fig.[4ji, exhibited this behavior. 

Two limiting cases of zero friction and complete friction at the specimen-platen 
interface were modeled. The simulation was performed by pushing the rigid platen 
increasingly to displacement of 0.485 mm, which is e = 0.02 deformation for the zero 
friction case. The effect of extensometer gauge length was tested by following the 
axial displacement of two nodes. Figure [l3j located l = 6.35 mm and l = 12.7 mm 
above the center plane of the specimen. These nodes correspond to the contact points 
of extensometers with G = 12.7 mm and G = 25.4 mm respectively. To calculate 
the stress-strain behavior, the force on the platen, the displacement of the platen, the 
displacement of the two extensometer contact points and a surface corresponding to the 
length between the platens, see Figure[l3] were recorded. The stresses were calculated 
from the measured radii at those points. 

Figure [14] shows the computed true stress-strain curves for strain measures across 
three gauge lengths (G = 12.7 mm, G = 25.4 mm, and G = 48.5 mm-across the 
loading platens) and for the two limiting friction cases: frictionless (/ = 0) and perfect 
friction (/ = 1). It also shows the true stress-strain curve used as the input model. In 
the figure, the true stresses, are calculated from the load P and the instantaneous area, 
A, by by assuming conservation of volume in deformation: 


AqIq = Al 


(9) 


where A (l and l 0 are the original area and gauge length respectively. This calculation 
mimics what would be done in analysis of an actual experiment, where the instanta- 
neous load-bearing area, A is calculated from the strain. For the frictionless case, all 
the strain measures yield the same stress-strain behavior as the input model, as they 
should. For the case of perfect friction, / = 1, the strains measured across the two ex- 
tensometer gauge length, G = 12.7 mm and G = 25.4 mm differ by less than 0.5 MPa 
each other and from the frictionless cases. The difference between the average of 
yield strengths measured by laboratories that lubricated and those that did not was was 
twenty times this large: A 75(0.2 % offset)=10.4 MPa, Figures [IT] and 12 The true 
stress calculated from the displacement of the loading platen, a common method that 
use the displacement of the actuator, corrected for machine compliance, is about 6 MPa, 
or 1.4 %, above all the other curves. Note that no laboratory measured strain this way; 
each laboratory used an extensometer on the test specimen, with contact points remote 
from the test specimen end faces. The agreement between the extensometer methods 
for the frictionless and complete friction cases eliminates friction as the source of the 
difference in yield behavior between the laboratories that lubricated and those that did 
not. 


The finite-element analysis also determined the evolution of the stress distributions 
in the test specimen. Figure [15] shows the von Mises stress evolution for four strains. 
Stresses less than the yield stress of the material model, a = 190.4 MPa, are shown 
in white. Consistent with other analyses |29ll6l, a non-deforming cap forms at the test 
specimen end. As the specimen deforms, a region of slightly higher stress develops 
that extends just into the region that the longer, G = 25.4 mm extensometer samples. 
Figure [16] shows the evolution of the shear stress, cri 2 for the same four strains. The 
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bulk of the test specimen has zero shear stress. Only at the test specimen ends, where 
friction constrains the deformation, does shear strain develop. 

4.3 Elastic Modulus 

The elastic modulus measured in the test is very sensitive to the alignment and sensitiv- 
ity of the test setup. As a result, ASTM E 9 12] requires that laboratories validate their 
compression test fixture by measuring the elastic modulus of aluminum alloy 2024- 
T4 bar, according ASTM E 1 1 1 ED and obtaining measurements on five consecutive 
specimens to within 5 % of a stated value, E= 73.8 GPa. The behavior of the reported 
elastic modulus should be a good indicator of the quality of the compression test facil- 
ity. 

The interlaboratory study did not require laboratories to revalidate their test setups, 
but they did report the measurements of elastic modulus, Figure [5] made during the 
compression tests of this study. Most of the laboratories only used a ASTM E 83 |32l 
Class B-2 extensometer, see Table [3] rather than the roughly twice-as-accurate class 
B-l extensometer required by ASTM E 1 1 1 eh when elastic modulus is the primary 
parameter to be determined. But because the primary parameter to be determined was 
the yield strength, E 9 Section 9.2 permits the use of the class B-2 extensometer. 

The laboratories had wide latitude in the method to determine the elastic modulus 
to report. To reduce the variability due to the interpretation of the individual labora- 
tory, and to put each modulus measurement on a consistent basis, we reanalyzed each 
stress-strain curve to estimate the elastic modulus using a technique ll33l based on the 
method of Scibetta and Schuurmans J34). These reanalyzed moduli are denoted by the 
symbol Eqa , to differentiate them from the moduli that the participants calculated and 
reported, denoted by E. In addition to estimating the modulus, the method includes 
some checks on the quality of the data that may be useful for laboratories that seek to 
improve their modulus measurement. 

Many of the following plots evaluate the quality of the measurements by employing 
the absolute modulus error, AEqa* in a measurement, which is defined as the absolute 
value of the difference between the accepted modulus l25l for 2024-T351 aluminum, 
Eai = 7b .‘2 GPa and the calculated modulus, TV ; .4 ■ 

AE G a = \E^ C - E ga \ GO) 

Other references [ 35] j36] [37] [2 ] recommend different values in the range 73.8 GPa < 
E < 76.1 GPa, and differentiate lf36l between longitudinal and transverse orientation. 

Method for analysis The method consists of nine steps. 

1 Find the knee of the stress-strain ( s — e) curve by the method of Scibetta and 
Schuurmans |[34| . 

1. Find the point on the stress-strain curve (ei, si) closest to si = 0.05max(s) 

2. Create point (e 0 , s 0 ) where e 0 = ei and s 0 = 0.2max(s). 

3. Find the point on the stress-strain curve ( et,St , ) that is on a line drawn 
through (eo, sq) that is tangent to the stress-strain curve. 
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2 Truncate the stress-strain data at (et, St) and normalize by these values. Retain 
the maximum values of the test record: (e max , s max )- 

3 Check digital resolution, S, of stress and strain. The digital resolution should be 

S e < 2?2 and 5 S < No more than more than 25 % of points should 

have zero stress or strain change from the previous point. 

4 Use the Scibetta-Schuurmans (34] algorithm to determine optimum fit window 
and calculate the slope, which is the elastic modulus -Ec^in normalized form. 
The Scibetta-Schuurmans algorithm finds the optimal region for determining the 
modulus by performing a linear regression of stress on strain on every possible 
subset of data in the region e 0 < e < e< that contains at least 20 % of the data. 
The regression with the lowest residual standard deviation is the optimum fit. 
In some cases, more than one million fits were evaluated for each stress-strain 
curve. 

5 Check for excessive noise in optimum fit window, defined as the standard error 
of the regression in step 4: Sj defined in Table[6]where n is the number of stress- 
strain pairs in the optimal region, Y, 3 is the normalized stress evaluated at point 
j, and Yj is the predicted normalized stress from the fit to the optimal region. 
Perform the complementary regression of strain upon stress as well. In both 
cases, the standard error in strain or stress should not be greater than 0.01. 

6 Extend the range of fit to include all stress points whose deviation from the opti- 
mal fit line is less than one times standard error computed in Step 5. 

7 Refit the extended data set to determine modulus, Eqa- 

8 Examine the shape of the stress and strain residuals as a function of strain in 
the extended fit. If the slope of the residuals in the first or fourth quartile of the 
extended range is more than 0.05, the data exhibit excessive curvature. At least 
five points must exist in each quartile to evaluate this curvature. 

9 Check that the stress range of the optimal fit is greater than 0.4st. 

The method employs four metrics to determine the quality of the modulus: 

• Data quality 1: digital stress and strain resolution should be sufficient (Step 2) 

• Data quality 2: strain or stress signal should not have excessive noise (Step 5) 

• Fit quality 1: the stress and strain residuals in extended data range should not 
have excessive curvature (Step 8) 

• Fit quality 2: the size optimal fit range, Rf, should be greater than or equal to 
40 % of the total range used to determine the modulus (Step 9) 

The rest of this section will examine the results of the modulus measurement, Eqa, 
using these four quality metrics. 
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Data Quality 1: sufficient digital resolution Figure [17] summarizes the strain and 
stress resolution. Estimating the digital resolution of a data set post-test is difficult, and 
the interested reader should consult the original reference by Graham and Adler ll33l . 
To create the figure, the digital resolution was estimated from the data set by taking the 
absolute value of the difference of each sequential pair of stress and strain points, for 
example |sj — Sj+i|. Those differences, expressed in units of 6 from step 3, were binned 
to create a histogram. If too many of the points show no change between reading, or if 
the most frequent change is many times larger than the digital resolution, S defined in 
Step 3, the stress-strain curve may have insufficient digital resolution. This condition 
could arise, for example, if a load cell were used near the very lower limit of its working 
range, or if the analog-to-digital converter had insufficient precision. Figure [TT] plots 
the index of the bin with the maximum fraction of the data, in units of 6, defined in 
step 3, against the fraction of the points in the “zeroth” bin, which is the fraction of 
points where the stress or strain value did not change between readings. The symbols 
show the fraction of the points in the bin of maximum fraction, broken into two groups: 
those where the bin of maximum fraction contained less than 25 % of the data ((0, 25]) 
and those that contained more than 25 % of the data ((25, Inf). The method identifies 
experiments with insufficient digital resolution as those where the bin of maximum 
fraction is greater than 3 and either the fraction in the zeroth bin is > 25 % or the 
fraction in the bin of maximum fraction is > 25 %. The first condition is identified by 
the dashed lines. The second “or” condition comprise those points identified by green 
circles with y value greater than 3. No experiment demonstrated insufficient digital 
resolution. 


Data Quality 2: excessive noise Excessive noise in the strain or stress signal will 
degrade the quality of the measured modulus. The method identifies experiments with 
excessive noise in these signals as ones in which the standard error of the fit is greater 
than 0.01. Figure p~8] plots the standard error of the fit for both strain and stress by 
laboratory. The symbols identify three levels of modulus error, Eq. [10] The noise in all 
the experiments was sufficiently low. No obvious relation exists between the absolute 
modulus error, A E, and either the standard error of stress or strain. 


Fit Quality 1: Excessive curvature in extended data set Some metric of the curva- 
ture of the optimal fit range is necessary, since the method only selects the best fit over 
an optimal region. That fit might still be poor. One method for examining the quality 
of the fit is to examine the deviation from the fit line in the first and fourth quartiles of 
the optimal range. Curvature of the stress-strain record frequently appears this way. If 
the slope of the residuals vs. strain in the outer quartiles is greater than 5 % of the slope 
in the optimal region, the curvature of the fit is deemed excessive. 

Figure[l9]shows the slope of the residuals in the outer quartiles. Different symbols 
show the level of absolute error of the elastic modulus, Eq. 10 All but one of the tests 
fall inside the ±5 % limits. In addition, linear regression of the absolute modulus error 
against the absolute value of the residual slope reveals no relationship. 
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Fit Quality 2: Final fit range The second measure of the quality of the fit to deter- 
mine the elastic modulus. Eqa, is the size of the optimal fit range, Rf. Figure |20|plots 
the absolute modulus error, Eq. |T0} against the fraction of the range of original stress 
data. The dashed line is the boundary of the minimum acceptable limit Rf > 0.4. 
The absolute modulus error, Eq. |T0} increases as the size of the fit region decreases, as 
expected, and the boundary between the two regions lies at the chosen final fit range 
minimum Rf = 0.4. The mean error for tests with fit region Rf < 0.4 is 2.5 times 
larger (4.88 GPa vs. 1.95 GPa) than that from the acceptable region. 

Conclusions from elastic modulus analysis Of the quality metrics, the size of the 
final fit range is the best predictor of overall quality of the modulus measurement, and 
the level chosen, Rf > 0.4, is a good metric for identifying potential problems with 
modulus measurement. 


5 Conclusions 

Five conclusions can be drawn from this study. 

• The repeatability of yield strength determined from compression tests conducted 
according to ASTM E9 can be expected to be about 1.1 % of the mean value, 
cv r = 0.01 1 , Figure[6]and Table[5] 

• The reproducibility of yield strength determined from compression tests con- 
ducted according to ASTM E9 can be expected to be about 2 % of the mean 
value, cvr= 0.020 , Figure[6]and Table[5] 

• Despite the perceived difficulties with alignment in compression, the repeatabil- 
ity and reproducibility of the compression test was among the best measured for 
uniaxial tests, Figure[9] 

• Using two diametrally opposed extensometers instead of a single extensometer 
can improve the precision of the strain measurement by more than two times. 
Figure [TO] 

• If the final fit range for estimating the modulus is less than 40 % of the elastic 
region, the modulus measurement is frequently in error. Figure [20] 

In addition, although the reported data indicated that yield strengths measured on unlu- 
bricated conditions were 3.1 % lower than those measured using lubricated specimens. 
Fig. [12] a finite-element analysis of the deformation that incorporated test-specimen- 
end friction did not identify any mechanism by which friction could caused the dis- 
crepancy. 
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Figure 10: Reported (a) modulus and (b) 0.2 % offset yield strength, YS (0.2 % offset), 
identified by number of extensometers used. 
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Figure 11: Effect of lubricant on stress-strain behavior. 
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Figure 12: Effect of lubricant on (a) modulus E and (b) yield strength 75(0.2 % offset). 
Dashed lines are the mean values for all data, including “not-reported,” not lubricated 
(“No”) and lubricated (“Yes”) tests. 
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Figure 13: Finite element mesh and boundary conditions of the upper-right quadrant 
of the test specimen. 
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Figure 14: Computed true stress-strain curves for the cases of frictionless (/ = 0) and 
perfect friction (/ = 1) and different gauge lengths, G: (a) full curve, and (b) expanded 
section that shows the small differences between the results at large strain. 
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Figure 15: Von Mises stress distribution in the upper-right quadrant of the test specimen 
for four true strains: (a) e = 0.002, (b) e = 0.008, (c) e = 0.015, and (d) e = 0.022, 
measured from the G = 12.7 mm extensometer for the case of perfect friction. 



Figure 16: Shear stress distribution in the upper-right quadrant of the test specimen for 
four true strains: (a) e = 0.002, (b) e = 0.008, (c) e = 0.015, and (d) e = 0.022, 
measured from the G = 12.7 mm extensometer for the case of perfect friction. 
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Figure 17: (a) Strain and (b) stress data resolution. Dashed lines enclose the region of 
insufficient resolution. Symbols denote the fraction of points in the bin of maximum 
fraction. 
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Figure 18: Noise in calculated modulus based on (a) strain, and (b) stress, calculated 
in Step 5. Symbols denote three increasing levels of absolute modulus error, Eq. [T0| 
All measurements are below the acceptance limit of 0.01. Data have been jittered to 
prevent over-plotting of points. 
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Figure 19: Slope of residuals in the (a) 1st, and (b) 4th quartile of the extended fit for 
modulus. Symbols denote three increasing levels of absolute modulus error, Eq. [10] 
Dashed lines denote the acceptance limit Data have been jittered to prevent over- 
plotting of points. 
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Final fit range, fraction 


Figure 20: Absolute error in modulus, Eq. 10 as a function of the fraction of the final fit 
range, Rj, used to calculate the modulus. Symbols identify the different laboratories. 
Points to the left of the dashed line are deemed to have insufficient fit range. 


